Backaction of a driven nonlinear resonator on a superconducting qubit 
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We study the backaction of a driven nonlinear resonator on a multi-level superconducting qubit. 
Using unitary transformations on the multi-level Jaynes-Cummings Hamiltonian and quantum optics 
master equation, we derive an analytical model that goes beyond linear response theory. Within the 
limits of validity of the model, we obtain quantitative agreement with experimental and numerical 
data, both in the bifurcation and in the parametric amplification regimes of the nonlinear resonator. 
We show in particular that the measurement-induced dephasing rate of the qubit can be rather small 
at high drive power. This is in contrast to measurement with a linear resonator where this rate 
increases with the drive power. Finally, we show that, for typical parameters of circuit quantum 
electrodynamics, correctly describing measurement-induced dephasing requires a model going beyond 
linear response theory, such as the one presented here. 

PACS numbers: 85.25.Cp, 74.78. Na, 03.67.Lx, 42.50.Lc, 42.65.Wi 



I. INTRODUCTION 

Two-level systems (TLS) and harmonic oscillators are 
the two simplest systems that can be described exactly 
with quantum mechanics. Consequently, many physical 
systems are described at least approximately by either 
of these two building blocks. As an example, in cavity 
quantum electrodynamics (CQED) pQ, an atom, mod- 
eled as a TLS, interacts with a photon field inside a high 
quality optical or microwave resonator, modeled as a har- 
monic oscillator. Another example is circuit quantum 
electrodynamics (cQED) [2 , cavity QED's little brother 
and a promising candidate for the realization of a future 
quantum computer [3]. In circuit QED, a superconduct- 
ing artificial atom (or qubit) [4] is coupled to a coplanar 
waveguide resonator. In the context of quantum informa- 
tion processing, the resonator both acts as a filter, partly 
protecting the qubit from decoherence and relaxation, 
and as a measurement device for the qubit state. 

However, contrary to cavity QED where the atomic 
properties are fixed, the engineered devices studied in 
circuit QED can be tuned and are custom built. There- 
fore, while devices dating from the early stages of circuit 
QED [2] [5] were well described by two-level systems cou- 
pled to harmonic oscillators, more recent qubits, such as 
the transmon [BHH] , the low impedance flux qubit [9] , and 
the tunable coupling qubit [TO] are better described by 
multi- level systems (MLS). This is also the case for the 
phase qubit [TT]. Moreover, while the standard architec- 
ture for qubit readout has long been linear resonators [2], 
many recent results |12H16j now use resonators made non- 
linear with embedded Josephson junctions. Not only do 
these nonlinear resonators provide a bifurcation amplifier 
regime which considerably improves the readout — a key 
requirement for quantum information processing — but 



they also exhibit remarkably enriched physics. As exam- 
ples, they have been used to parametrically amplify small 
signals [HI [18] and generate squeezed light [19] . 

The performance of nonlinear resonators as parametric 
amplifiers for small signals [20] as well as their backaction 
on a qubit have also been studied theoretically [5T] [55] . 
However, in Refs. [HI [55], the qubit was assumed to be 
a two-level system, something which is often insufficient 
to understand many types of superconducting qubits. 
Moreover, a linear response of the output signal to the 
input (qubit) signal was assumed. While linear response 
holds away from the nonlinear resonator's critical point, 
where bifurcation becomes possible, and away from the 
switching thresholds in the bifurcation amplifier regime, 
we show that it breaks down close to these points. We 
show that linear response is unlikely to be sufficient to 
describe a qubit readout with a nonlinear resonator when 
considering typical cQED parameters. Finally, the usual 
dispersive theory with linear resonators assumes driving of 
the resonator close to its resonance frequency for measure- 
ment [2, 23 . As a result, the theory obtains a dependence 
of the ac-Stark shift on the frequency detuning between 
the qubit and the resonator, rather than between the 
qubit and the measurement drive. This is especially im- 
portant when measuring with a nonlinear resonator since 
there is always a significant frequency detuning between 
the drive and the resonator in such cases. 

In this paper, we derive a reduced qubit model go- 
ing beyond these assumptions. We do so using unitary 
transformations, especially the dispersive 23, 24 and the 
polaron transformations [55H5T] . We are especially inter- 
ested in describing the ac-Stark and Lamb shifts of the 
qubit as well as its measurement-induced dephasing |28j . 
We note that this theory was developed in parallel to 
and already tested against the experimental results of 
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FIG. 1. (Color online) Representation of one possible im- 
plementation of the system considered in this paper. This 
represents a stripline resonator (blue) made nonlinear with an 
embedded Josephson junction (dark green), capacitively cou- 
pled to a transmon qubit between the central conductor and 
the ground planes. The model described in this paper however 
applies to various other nonlinear resonators and qubits (see 
text). 

Ref. [H]. 

In section [TTJ we write the general master equation that 
is used to describe the multi-level qubit coupled to the 
nonlinear resonator. In section |ITl| we recall the minimal 
multi-level system model of linear circuit QED in the 
dispersive regime. In section [TV] we describe the basic 
characteristics of nonlinear resonators and explain why we 
need to go beyond the assumptions given in the previous 
paragraphs. In section [Vj we derive a reduced model for 
the qubit through a series of unitary transformations. In 
section |VI| we compare the predictions of the analytical 
model to experimental |16] and numerical data and find 
quantitative agreement within the limits of the model. We 
also explain how the ac-Stark and Lamb shifts as well as 
the measurement-induced dephasing are changed by the 
nonlinearity of the resonator. We finally test the regime 
of validity of the linear response theory and show that 
it is unlikely to be sufficient to describe any high-fidelity 
qubit readout with a nonlinear resonator. 



II. PRESENTATION OF THE SYSTEM 

We consider a system made of a multi-level qubit cou- 
pled to a nonlinear resonator. We describe the nonlinear 
resonator with the Hamiltonian (h — 1) [5D] 

4- I£ 4. 4. a.3 rj , 

H r = uj r a ] a + —a a aa H — — a 1 a , (2-1) 

£ o 

where a^' are the annihilation (creation) operators, uj t 
is the resonator low-power resonance frequency, and K 
and K' are quadratic and cubic Kerr constants. Such a 
Kerr nonlinear resonator could be an LC-circuit with an 
added Josephson junction |12j or a stripline resonator with 
one [H] (see Fig. [1} or many [T71 [55] embedded Josephson 
junctions. In all these cases, the Josephson junctions act 
as nonlinear dissipationless inductances, rendering the 
resonator nonlinear. 

We describe the qubit by the generic many-level system 



Hamiltonian 

M-l 

H q =Y. ^n li4 = n w , (2.2) 

i=0 

where M is the number of qubit levels, u>i is the frequency 
of the qubit eigenstate Ujj = \i) and where we 
have introduced the short-handed notation 

JW-l 

n x = x i u i,h ( 2 - 3 ) 

i=0 

which we will use on multiple occasions throughout this 
paper. The eigenstates {\i}} could be for example charges 
tunneling on and off a superconducting island such as 
for a Cooper-pair box [3D], superposition of such charges 
for a transmon qubit (SJ or current flowing clockwise or 
counterclockwise in a superconducting loop for a flux 
qubit |5T] . 

We assume a dipolar coupling between the qubit and the 
resonator and describe it by the interaction Hamiltonian 

M-2 

Hl =J2 #( ot + «)( n M+i + Hi+i.i), (2.4) 

i=0 

where gi are the coupling constants. The only constraint 
on the qubit that we impose for our model is that the 
selection rules only allow transitions between the qubit 
states \i) and \i ± 1) through the resonator. This restric- 
tion is fullfilled for good two-level qubits such as the 
Cooper-pair box [30], the phase [32] and flux [31] qubits, 
but is also realized for some more recent multi-level qubits 
such as the transmon [DJ [33J and the low impedance 
flux qubit 0. 

To understand the experiment of Ref. [16] . we also 
consider driving of the resonator. We allow for multi- 
ple qubit-detuned drives d £ {d±, d,2, d n } as well as 
one spectroscopy drive s, quasi- resonant with the qubit 
frequency, that we model by the Hamiltonians 

H d = W^ct + eje^o, (2.5a) 

d 

H s = e^^a) + e*e^*a, (2.5b) 

where and €d, s and LOd.s are the drives' amplitude and 
frequency. By quasi-resonant, we mean that uj s is always 
much closer to the |0) <H- |1) qubit frequency than to any 
other qubit transition frequencies. In experiments, these 
drives take the form of microwave signals sent to one port 
of the resonator and either transmitted to the other port 
or reflected back depending on the circuit design. As 
in the experiment of Ref. [16] , we will later on take the 
amplitude of the spectroscopy drive e s to be small such 
that its contribution to the intra-resonator field is small. 
The case of high amplitude spectroscopy will be treated 
in a following publication [53) . 
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Finally, to model dissipation, we use the Lindblad-type 
master equation 

p = —i [H, p] + nT>[a]p + [a 2 ]p 

M-2 / \ 2 (2 6) 

+ 7 E - t> [n M+1 ] P + 2 lv v [n £ ] P , 



where 



-(2ApA^ - A^Ap - pA^A), (2.7) 



. In this master 



/ d(cji-^o) ~\ 1 



and = H r + H q + Hi + H d 
equation, n and knl are the resonator's rates of one- and 
two-photon loss [2D], 7 is the qubit |1) — > |0) decay rate 
and 7^ is the qubit pure dephasing rate for the same 

states. For Il e , we defined £i = d(^-"o) 

as the X-dispersion, where X is some control parameter 
(could be flux or charge for example), with £0 — and 
£\ = 1 by definition. 

This master equation can be obtained by modeling 
the coupling of the qubit and the resonator to baths of 
harmonic oscillators and then tracing over the baths [34] . 
When obtaining this master equation, we made three 
assumptions. First, we assumed that the noise spectra 
are white around the relevant frequencies for relaxation 
(~ GHz) and dephasing (< 1 MHz). For this approxima- 
tion to hold, the baths must be white on a frequency range 
comparable to the resonator or qubit linewidths. While 
this approximation should hold for relaxation (~ GHz fre- 
quencies) if the resonator and the qubit have high quality 
factors, it may fail for dephasing (< 1 MHz frequencies) 
if, for example, the noise has a 1// spectrum and hence 
varies by many orders of magnitudes over a single res- 
onator or qubit linewidth. In this latter case, one needs to 
be more careful and take the noise spectrum into account 
when deriving the master equation [6, 35.. Second, we 
assumed that the noise causing qubit relaxation couples to 
the qubit through dipolar interaction, yielding the scaling 
in gi/go for the 7 dissipator. Finally, we considered that 
dephasing is caused by (white) noise at low frequencies 
in the control parameter X. 



In this regime, there is no direct exchange of energy be- 
tween the qubit and the resonator, and most of the physics 
can be understood from an approximate diagonalization 
of the undriven Hamiltonian H r + H q + Hi [2] . To second 
order in perturbation theory and assuming that the qubit 
is a TLS, this diagonalization yields 



-^disp — 



X , t , t 

— (7 z + u> r a ] a + x&za a, 



(3.1) 



where the effective qubit frequency is Lamb-shifted by 
a quantity \ — g 1 /A^o- The last term of this Hamil- 
tonian can either be seen as a qubit state-dependent 
pull of the resonator frequency — which allows for qubit 
measurement [2] — or as an ac-Stark shift of the qubit 
frequency that depends on the number of photons in the 
resonator [42]. 

In addition to the Lamb and ac-Stark shifts of the 
qubit frequency, the qubit 's coupling to the driven res- 
onator leads to additional sources of relaxation and de- 
phasing. Among these are Pureed relaxation [8] in which 
the qubit relaxes through the resonator's photon loss chan- 
nel, dressed dephasing [23] |33] H3] in which pure dephas- 
ing of the dressed qubit-resonator states leads to effective 
relaxation and heating of the qubit, and measurement- 
induced dephasing [27, 28, 42J which is the unavoidable 
dephasing caused by acquisition of information about a 
quantum system. For a linear resonator and in a disper- 
sive measurement regime, it is shown in Refs. [27, 28] that 
the measurement-induced dephasing rate is given by 



_ kD 2 k_ 



(3.2) 



where D = \a\ — ao| is the distinguishability of two pointer 
states of the resonator and n is the average number of 
photons inside the resonator. Under resonator driving, the 
pointer state at is the coherent state \ai) that represents 
the resonator's field if the qubit is in the state \i). For 
a linear resonator and a two-level system described by 
the dispersive Hamiltonian Eq. (3.1) with a single added 



drive of amplitude e p and frequency uj p , these coherent 
states are given by 



III. LINEAR CIRCUIT QED IN A NUTSHELL 

Before going to the nonlinear case, it is useful to review 
some aspects of the more standard linear case. In linear 
circuit QED, one is interested in the system described 
in section [TlJ but with K = K' = knl = 0, and with 
a qubit which can have two or more states. Many as- 
pects of this system have been studied extensively both 
theoretically and experimentally, ranging from qubit mea- 
surement 2 , 36 and single- and two-qubit gates [S] |3"TH4"T] 
to dissipation and dephasing [8] 123 ] [27 ] [28 ] I42H44] . In this 
section, we present the minimal theory of the dispersive 
regime where the couplings gi are much smaller than the 
qubit-resonator detunings A^., = Uij — uj r = lji — ujj — uj r . 



ai/o 



—i(uj r — uj p ± x) ~ K /2 ' 



(3.3) 



and are represented in phase space on Fig. [2] for a resonant 
drive. 

The distance D between these pointer states in phase 
space depends on the cavity pull x and, for a dispersive 
measurement with a linear resonator, increases with the 
number of photons or equivalently with the strength of the 
measurement drive. It is further shown in Ref. [27] that, 
in the linear case, the measurement-induced dephasing 
rate reaches the smallest value permitted by quantum 
mechanics. In other words, it saturates the inequality 



r > 



r,. 



(3.4) 
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FIG. 2. (Color online) Phase space representation of the 
pointer states ao (full blue circle) and ai (red dashed circle) 
given by Eq. (3.31. Parameters are x/2t = 4 MHz, k/2tv = 
10 MHz, lo p = Lo r and e p /27r = 30 MHz. 



where r meas . is the measurement rate [45 , corresponding 
to the rate at which information is gained on the system 
being measured. One of the questions that we will try 
to answer in this paper is whether or not this inequality 
can be saturated when using a nonlinear resonator for 
homodyne dispersive measurement of the qubit. 




FIG. 3. (Color online) (a) Amplitude of the resonator inter- 
nal field (arbitrary units) in response to a drive of reduced 
frequency fl = 2(uj r — bid) In (S^c = v3) for increasing drive 
amplitudes e^. The horizontal lines indicate regions for which 
linear response theory would (full green lines, check marks) and 
would not (dashed red lines, X marks) be valid for modeling a 
qubit-resonator system, (b) Stability diagram of the resonator 



(see text for explanations). As discussed in subsection VIA 
the vertical dashed lines represent the two operating points 
studied in this paper. 



IV. FEATURES SPECIFIC TO NONLINEAR 
CIRCUIT QED 

Depending on the amplitude Cd and frequency cud of 
the drive, the response of a Kerr nonlinear resonator can 
be either mono- or bi-valuated. The stability diagram 
describing this behavior can be parametrized by the re- 
duced detuning frequency fl = 2{tu r — LJd)/ K and by the 
drive amplitude e<j. If the reduced detuning is smaller 
than — but close to — a critical value flc = \/3, the 
nonlinear resonator can be used as a low-noise parametric 
amplifier |17j . This has been used recently to amplify 
microwave signals at the single photon level |18] , For 
0,/Q.C > 1, the stability diagram, illustrated in Fig. ^b), 
shows two bistability thresholds [46 . Below the first one 
(dashed green line), a low (L) amplitude response [see 
Fig. [3^a)] of the resonator is observed. Above the sec- 
ond one (full red line), one rather observes a high (H) 
amplitude response. Between the two thresholds, both L 
and H are stable. Because of the coupling to the qubit, 
this stability diagram depends on the qubit state. This 
dependence allows the nonlinear resonators to be used as 
a sample-and-hold detector as has been demonstrated in 

Refs. pa nam]. 

Before going forward with the theory, we want to high- 
light two peculiarities of circuit QED with a nonlinear 
resonator that are often overlooked. These two aspects 
— the detuning of the readout drive from the resonator 
frequency and the limits of the linear response theory — 
as well as their impact on the theory are discussed further 



in the following subsections. 



A. Detuned measurement drive 

Both in usual low power dispersive measurement of 
a TLS [2 j and the more recent high power avalanche 
readout [48-50], measurement with a linear resonator 
is done with a drive at or very close to the resonator 
frequency uj r . On the contrary, measurement with a 
nonlinear resonator is always done with a drive source 
significantly detuned from w r [131 [TSl S3 [5T] . As can be 
seen from in Fig. [3j this detuning is required to bias the 
system either in the region of highest parametric gain or 
in the bistability region. 

Because of the Jaynes-Cummings interaction, the drive 
on the resonator also acts on the qubit. Since the cavity 
is acting as a filter, the effective drive amplitude as seen 
by the qubit is expected to scale as l/(u> r — UJd)- Photons 
entering the cavity because of this drive will cause an 
ac-Stark shift of the qubit x( a ^ a )- The shift per photon 
X should depend on the drive frequency. This is however 
not the case for the usual expressions for a TLS, where 
X = Xo 0, or for a MLS, where \ = 2Xo~Xi M> witn Xi = 
gf/A i+ i^i — gf/(uj i+ i :i — io r ). Indeed, these expressions 
scale with the inverse of the qubit-resonator detuning. 
One would rather expect to find \i = ffiV( w j+M — u <i) 
since the drive photons are at frequency UJd- While a 
relative change of a few percents on Ai.o yields the same 
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relative change on x for a the two-level system, the effect 
can be twice as big for a MLS because of the reduced 
value of x- To obtain quantitative agreement with the 
results of Ref. [TB] , we obtain below an expression for the 
ac-Stark that contains the expected frequencies. 



B. Limits to the validity of linear response in 
circuit QED 

As stated before, Kerr oscillators have been used ex- 
perimentally as parametric amplifiers for small signals. 
They have also been studied theoretically extensively. As 
examples, Yurke and Buks have studied their performance 
and calculated their gain [50] , while Laflamme and Clerk 
have shown that these amplifiers are not quantum limited 



in the sense of Eq. (3.4 1 for a qubit measurement 



Moreover, these last authors show that the quantum limit 
can be reached if one makes use of correlations between 
the resonator and the system coupled to it. 

These two results were however obtained in the limit of 
linear response theory. In this limit, one finds the driven 
resonator's stationary state a without the coupling to the 
qubit and then expands the solution including the qubit 
around the stationary solution a ~ a + 8a. For a qubit 
measurement, the signal that is amplified by the resonator 
takes the form of a pull ±x °f the resonator frequency 
which in turns depend on the qubit state as expressed in 
Eq. (3.1). For a linear resonator in the dispersive regime, 



the ccj's given by Eq. (3.3) can be rewritten as 



«i/o = a I 1 
wall 



—i(uj r — tjjp ± x) ~ K /2 
±iX 



—i(to r — ujp) — n/2 J ' 



(4.1) 



where the linear response expressed by the second line 
holds if \i(cj r — oj p ) + k/2\ ^> \x\- Therefore, the validity of 
linear response in this linear dispersive case is not affected 
by the driving strength, but is rather determined by the 
ratio 2x/k <C 1. 

This analysis however does not hold for a nonlinear 
resonator. Indeed, in order for linear response theory 
to stay valid with a nonlinear resonator, a must change 
linearly with the pulled frequency — or equivalently with 
the drive-resonator detuning — over a frequency range 2%. 
While for a linear resonator, it has been shown that 
the optimal SNR is obtained for 2x — K, the improved 
measurement efficiency with a nonlinear resonator allows 
for smaller cavity pulls. Taking x = 0.2k as a typical value 
of the cavity pull translates into a range of fl/Qc ~ 0.5 
over which the signal must be linear in frequency for the 
linear response to stay valid. This range is illustrated 
on Fig. [3](a) with the horizontal lines. The full green 
lines represent regimes for which linear response would be 
a good approximation, while dashed red lines represent 
regimes for which the response is not linear over the 
appropriate range. We argue that the linear response 



approximation will break down as soon as the slope of 
the response — and hence the gain of the amplifier — 
becomes significant. 

In the following section, we derive a theory that goes 
further than linear response theory using the polaron 
transformation approach of Ref. [2Z1 ■ 



V. REDUCED QUBIT MODEL 

In this section, we derive a reduced qubit model that 
captures the ac-Stark and Lamb shift of the qubit transi- 
tion frequencies as well as measurement-induced dephas- 
ing. This is done by performing unitary transformations 
on the master equation (2.6). These transformations have 
two objectives. First, transforming the system into its 



eigenbasis in which the physics is easier to understand. 
Second, eliminating the resonator to obtain a master 
equation for the qubit alone. 

In order to reach these objectives, many transforma- 
tions have been used in the litterature. The dispersive 
transformation [2H S3] (here generalized for a MLS) 



D 



exp 



M-2 

E 



(5.1) 



where A, = gi/{^i+i.i — oj r ), diagonalizes the Jaynes- 
Cummings Hamiltonian and reveals the Lamb and ac- 
Stark shifts. This transformation however only knows 
about photons that are at the resonator frequency ui r and 
fails to correctly model the measurement drive-resonator 
frequency detuning as discussed in Sec. |IV A| 

Another useful transformation is the displacement op- 
erator 152 , 



D(a) = exp [aa' — a* a] 



(5.2) 



which displaces a coherent state |— a) of a resonator to 
the ground state |0). In operator representation, it cor- 
responds to the change a — > a + a, where a represents 
the classical average field and a its quantum fluctuations. 
Doing this transformation before the dispersive transfor- 
mation, as was done for example in Ref. [38], yields the 
correct qubit-drive detuning in the ac-Stark shift. The 
ac-Stark shift then depends on the mean field amplitude 
a and the ac-Stark shift per photon depends on the drive- 
qubit frequency. However, doing this transformation in 
the context of a nonlinear resonator is akin to doing a 
linear response theory. Indeed, it is the same as assuming 
that the intra- resonator field is \a) and then look at all 
further perturbation, such as the cavity-pull, with respect 
to this mean field value. This will be discussed further in 

sec. inn 

A third transformation that was used in Ref. [27] to cal- 
culate the measurement-induced dephasing rate, as well as 
in Ref. 26 to study a qubit coupled to a mechanical res- 
onator beyond the rotating wave approximation (RWA), 
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is the polaron transformation [55] (here generalized for a find 
MLS) 



M-l 

i=0 



;D{ai). 



(5.3) 



This corresponds to a displacement transformation that 
is conditional on the qubit state. It allows for different 
cavity states jaj for each qubit state which makes 
it possible to go beyond the linear response approxima- 
tion. It is important to note that the field amplitudes m 
are free parameters in this transformation. In practice, 
these amplitudes will be chosen such as to cancel specific 
terms in the transformed Hamiltonian. Moreover, and 
as will become clear below, these different a^s will be 
independent solutions of qubit-state-dependent nonlinear 
equations, and not expansions around a mean solution a 
of a single mean nonlinear equation. 

In the following subsections, we perform three transfor- 
mations in order to approximately diagonalize the Hamil- 
tonian and transform the full master equation Eq. (2.6) 



to a reduced qubit master equation containing all the 
relevant physics needed to account the low power spec- 
troscopy of a qubit coupled to a nonlinear resonator driven 
by an external field. 



M-2 



n„ + £ ft [n^n M+1 + n I+M nj 

K 



i=0 
2 



im[n Q n* 



K' 



in r 



£ fce-^'n* +h.c], 



(5.5) 



where II Q is defined according to Eq. (2.3 1. In this Hamil 



tonian, the second term of the first line acts as drives on 
the qubit at the frequencies contained in the time depen- 
dence of a. The last two lines will be partly cancelled 
below by the choice of a given in Eq. (5.19) and we will 
neglect the small remaining parts. 

We also obtain the qubit-resonator Hamiltonian, limited 
to terms with one resonator ladder operator, 

M-2 

H[ = G'^+G'W J2 S i (a t +a)(n vi +i+IL, +M ), (5.6) 



i=0 



where 



C = K - if )n Q + (k - m NL )|n a | 2 n 



A"|n Q | 4 n Q + ^ 



— iuidt 



iU a . 



(5.7) 



A. Polaron frame 

While the polaron transformation can be performed 
exactly on terms that are diagonal in the qubit subspace, 
applying it on non-diagonal terms unfortunately yields 
complicated expressions. For example, applying it on a 
qubit ladder operator cr_ yields 



P = a—D(ai — aa)e 



-ilmfa? an 



(5.4) 



which, through the displacement operator D, contains 
all powers of a and aL For this reason, the polaron 
transformation was used in Refs. [23l [27l H3] after doing 
the dispersive transformation which eliminates the off- 
diagonal qubit operators. In this paper, we instead apply 
it before the dispersive transformation, assume that \a\ — 
cto\ <C 1 and take as a simplification PV P w er_. The 
small distinguishability approximation |ai — ao| < 1 will 
be made throughout this calculation and will limit the 
range of validity of the theory in a way which will be 
discussed later. 

The application of the polaron transformation on the 
master equation (2.6) is presented in Appendix |A| Follow- 
ing this Appendix, we use the notation H[ to represent a 
part of the Hamiltonian in this first transformed frame 
that contains i resonator ladder operators First, for 
« = corresponding to the qubit-only Hamiltonian we 



We will see below that the two first terms of H[ can be 
cancelled by a proper choice of n a . The last term will 
yield the Lamb shift of the qubit frequencies once the 
dispersive transformation is done. 

Finally, we find for the Hamiltonian containing terms 
with two resonator ladder operators 



where 



H' 2 = u' r (a)(Ja + (Ta t2 + T*a 2 ) 



/ r (a) = uj r + 2K\U a \ 2 + 3K'\U a \ 



^'|n Q | 2 )n 2 , 



(5.8) 



(5.9a) 
(5.9b) 



With |n Q | 2 = X^fio 1 l a i| 2 ni,i corresponding to the num- 
ber of photons associated to the different qubit states, 
we see from the expression for ui' r (a) that the resonator 
frequency is changed by the nonlinearity as expected. 
Moreover, the last term of Eq. (5.8) will squeeze the res- 
onator field. This will be studied in elsewhere [S3] and, 
for the scope of this paper, we will consider squeezing to 
be negligible. 

Having transformed the Hamiltonian, we now apply 
the polaron transformation to the dissipative parts of 
the master equation Eq. (2.6). We note that one could 



alternatively apply the transformations on the system- 
bath Hamiltonians before deriving the master equation. 
In this way, it would be possible to relax the white noise 
approximation |23j . something we will not focus on here. 
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Applying the transformation, we arrive at the master 
equation of the system in the polaron frame 



P 



-i [H' Q +H[ + H' 2 ,p'] + nV[a]p' + n Nh V[a 2 } P ' 

M-2 , y 2 

+ KV[u a ] P ' + K NL v[n 2 a } P ' + 4 KNL p[ a n Q ]p'. 

(5.10) 

When obtaining the dissipative terms, we have neglected 
non-Linbladian terms of the form a[//,n*] under the 
assumption that in the polaron frame, the resonator is 
in, or close to, its ground state (see Appendix [A] and 
Ref. [57] ). In this equation, the two first lines are the 
Hamiltonian part as well as the unchanged parts of the 
dissipative terms. The last line contains measurement- 
induced dephasing through the single-photon (first term) 
and two-photon (second term) loss decay channel, as well 
as some additional resonator decay (last term). 

In this polaron frame, we end up with a resonator whose 
frequency is shifted by the nonlinearity and the amplitude 
of the classical fields o^. This resonator is driven with 
an adjustable strength G which could be set to zero by a 
proper choice of a,-. It is important to note that we did 
not make that choice yet because, if we did, we would 
have cdj = oij and therefore would lose all dependence of 
the field amplitudes oti over the qubit state. The choice 
of the value of the qubit-state dependent fields an will 
be made only after moving, in the next subsection, to 
what we call the classical dispersive frame. Finally, in 
the polaron frame, the qubit is driven off-resonantly at 
frequencies uid and quasi-resonantly at frequency u s with 
amplitudes aj )( j and an s . As we will now show, the off- 
resonant drives will yield the correct ac-Stark shifts of 
the qubit frequencies. 



B. Classical dispersive frame 



We now focus on the qubit Hamiltonian Eq. ( |5.5[ ) Since where 



U a has a time dependence involving the drive frequencies, 
this Hamiltonian is that of a qubit driven with multiple 
direct drives. We have not yet computed the amplitude of 
the fields yet and we will do so now taking a = J2i a 



a quasi-resonant drive with amplitude giCti^ s and frequency 
lo s . Focussing for now on the drives w d ^ w s , the first line 
of this Hamiltonian can be approximately diagonalized 
with an analog of the dispersive transformation Eq. (5.1) 



exp 



M-2 

£ 

. i=0 



ern 



(5.12) 



where & is a classical analog of the operator AjCr . Because 
of this analogy, we will refer to this as the classical disper- 
sive transformation. This transformation is performed on 
the master equation ( |5.10[ ) in Appendix [B] where we take 



£6 



— iui d t 



(5.13) 



with £ itd = Afa it d. 

In the spirit of the dispersive transformation, Dc as- 
sumes an off-resonant driving and therefore cannot be 
applied to transform the spectroscopy drive s. When 
doing the transformation, we drop time-dependent terms 
involving two different drive frequencies ujd 1 ± ujd 2 under 
the rotating wave approximation. We also assume that 
for the purpose of getting the qubit transition frequencies, 
oti,d = Q!o,d- This is the same as taking |a, — oti+x\ to 
be small. Essentially, we assume that the difference in 
the pointer states is not important to describe the value 
of the qubit transition frequencies, but is important to 
describe their widths. In other words, we say that the 
mean transition frequency depends on the mean cavity 
field, which is approximately a$ d at low spectroscopy 
power (for the qubit is mostly in its ground state), while 
the width of the transition frequencies depend on the 
deviation of the cavity field from ao.d- 

Performing the above transformation on the qubit 
Hamiltonian H' to fourth order in perturbation theory 
together with the simplifications just outlined, we find 



M-l 

£ 

i=0 



M-2 

■ E ' 

i=0 



n, 



•h.c, (5.14) 



cjJ' = Wi + £;S?|a«,| 
d 



iE^i 



Ctd\ 



(5.15) 



-iui d t . 



This choice assumes that the are the ac-Stark shifted qubit frequencies and 



multiple drives are spread out enough in frequency such 
that one drive does not contribute significantly to the 
field oscillating at another drive's frequency. We therefore 
take 



M-2 



H' ^n u +J2J2 9i<Xi,de- iUdt n i+14 + h.c. 



if = -{xt-xU), 

If = -4S? (|Af | 2 + lAttl 2 ) 

- (3X? +1 \Aff - X?\Af +1 \ 2 ) 

+ (3Xt 2 \At 1 \ 2 -Xt 1 \At 2 \ 2 ), 



(5.16) 



i=0 d 



M-2 



(5.11) 



are the quadratic and quartic ac-Stark shift coefficients 
with 



i=0 



Transitions \i) O \i + 1) are then driven by an off-resonant 
drive with amplitude giOti id and frequency Ud, as well as by 



A" 



X' 1 



- u d 



(5.17) 



a>i+i — u>i — uJd 



We note that g { = for all i <£ [0, M - 2] in the ini- 
tial model such that terms with a negative index or an 
index above M — 2 on the right hand side of the equa- 
tions above vanish. Comparing these expressions with 
equations (3a) and (3b) of Ref. |49| , we highlight a few 
differences. First, both Sf and now depend on the 
drive frequency io d instead of the resonator frequency 
u r . As explained in section |IV A\ this follows from con- 
sidering that the driving photons can be at a frequency 
signif icantly detuned from v r . Actually, Eqs. (5.16) and 



(5.17) also hold for linear cQED, where the measurement 



drive is in practice chosen to be quasi-resonant with uj r . 
Next, the equation for Sf does not involve terms of higher 
order than g 2 . In Ref [49], these higher-order terms came 
from choosing a specific order for ladder operators when 
computing Ki (i.e. a^a^aa = a' aa^ a + a) . Here, the 
field is classical and there is no such ordering choice to 
be made. Finally, in Ref. [49], a second-order coupling 
caused by two-photon transitions was diagonalized, yield- 
ing fourth order corrections. This second-order coupling 
is however only significant in the straddling regime where 
the resonator frequency is between two qubit transition 
frequencies [B]. Since we are not considering this regime 
here, this two-photon transition was neglected. 

The next step is to apply the transformation Dp on 
H[ to find 



M-2 



H" « .9*( at + a)(n M +i + IT +M ) 



i=0 



(5.18) 



+ (G" a f + G' n a) + H 



SB, 



where G" can be found in Eq. ( B7 1 . The Hamiltonian H$b , 



whose definition can be found in Eq. (B8), corresponds to 
red and blue sideband transitions. This Hamiltonian is 
the multi-level equivalent of the one obtained in Eq. (BIO) 
of Ref. [33] for a two-level system driven by two detuned 
drives and experimentally studied in Ref. [54] . The drive 
strength G" can be set to zero with a proper choice of 
the fields , yielding an undriven resonator in this frame. 
Assuming that \cj dl — uj d2 \ is sufficiently large to neglect 
time-dependent cross terms, choosing G" = implies 



= {oJ r - u d - if )a ifd + (K - in-NL)\ui\ 2 ai,d 
+ K'\atfa lA + e d + (sf + ^Kf\ ai \ 2 ^j a lAl 

for each qubit-detuned drive and 

= (w r - lo s - if )a i s + (K - iK.^)\a.i\ 2 ai^ 



(5.19) 



K'\a,\ A a 



(5.20) 



for the spectroscopy drive. In writing these expressions, 
we have again assumed that, even though ctj ^ aij, these 
amplitudes are close enough to replace one by the other 
in order to uncouple the equations for i ^ j. 

We stress that because Eq. (5.19) contains the qubit- 



here go beyond the linear response theory for the response 
of the field to a change of the qubit state. As explained 
briefly in section [IV B| and as we will detail further later, 
a linear response theory would instead have solutions of 
the for m otj — a + /(Si), where a would be the solution 
of Eq. Km with Sf =Kf = and /(Si) would be some 



linear function of the cavity pull. 

We note that the equations for two different drives 
d\ 7^ e?2 are coupled through the total field |aj|. However, 
in the interest of reproducing the results of Ref. [16] , from 
this point on we will consider only a single qubit-detuned 
drive which we will label d = p (the pump drive) in 
addition to the spectroscopy drive s. This implies that 
only the first line of H'{ will remain. For the purpose of 
calculating a s , we will also assume that the spectroscopy 
amplitude e s is small enough so that \cti tP \ ^> |ai lS | and 



that cti, p ~ a.j 



such that we can replace «i by on 



in the equation for Oj )S . Finally, since in practice K 1 <C 
K <C 5, performing the classical dispersive transformation 
on H' 2 would yield corrections smaller than those that we 
have kept so far. We therefore neglect those and take 
H 2 = H' 2 . 

Finally, applying the transformation T>c on the dissipa- 
tion yields the master equation in this doubly transformed 
frame 



P 



-i [H'J, p"] + K "V[a]p" + 2 lv V TLJ p" 



M-2 



M-2 



J2 1UP [n M +i] p" + E P<+m] p" 

i=0 i=0 

kD [n Q ] P " + KNL v [ill] p" 

[M-l 



E 



90 



n, 



where 



(5.21) 

(5.22a) 

(5.22b) 
(5.22c) 

(5.22d) 



k" = k + 4K NL |a p | 2 , 
// _ ( 9i\ , ii 

7f,i = 7dd,z> 



and 



M-2 



state dependent cavity pull, the solutions a,- d obtained 



H'J = HS + E fliC^i+i + an i+ i,i) + Hg. (5.23) 



These two transformations result in an ac-Stark shifted 
qubit that is driven with a spectroscopy drive of amplitude 
a St i and frequency u„ coupled with a Jaynes-Cummings 
coupling to an undriven resonator whose frequency is 
shifted by the nonlinearity. This resonator sees addi- 
tional relaxation k" > k due to the two-photon-loss re- 
laxation channel. The qubit sees its intrinsic dephasing 
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at a rate 7^, as well as relaxation at rate 7"^ and heat- 
ing at rate 7"^ These relaxation and heating rates are 
modified by dressed-dephasing [33] (first term of 7dd,i)j 
but also dressed measurement-induced dephasing (sec- 
ond term). These rates were obtained assuming white 
noise for all the dissipation channels. If the noise is not 
white, the rate 7dd,i will depend on the noise spectra of 
the qubit dephasing and resonator relaxation channels at 
±(wj+i.i — ujd) [23 . In addition to intrinsic dephasing, the 
last two lines of Eq. (5.21 1 contain three other sources of 



dephasing. The first term will yield measurement-induced 
dephasing [25], while the second and the third represent 
respectively measurement-induced dephasing through the 
resonator two-photon loss decay channel and through the 
emission of an excitation by the qubit in its environment. 
While not measurable, this excitation in principle car- 
ries information about the qubit state and thus causes 
dephasing. 



C. Quantum dispersive frame and reduced master 
equation 

The final effect that we would like our model to capture 
is the Lamb shift of the qubit frequencies due to vacuum 
fluctuations of the resonator. To obtain this shift, we 



perform the dispersive transformation D of Eq. (5.1) 



on the master equation Eq. (5.21). Doing this while 



neglecting the photon population that is almost zero 
in the polaron frame [see discussion below Eq. (5.10)] 
yields the same master equation, but with the transformed 



Hamiltonian H' s " = H' " 



M—2 



H'i' with 



#o" = n*,'» + J2 ''<"<-' + h-c, (5.24a) 

i=0 

Hg' = [u)' r {a) + n s(Q) ] aU + To* 2 + TV (5.24b) 



and where the relaxation rate 



(5.25) 



has an added Purcell relaxation rate. The new Lamb- 
shifted frequencies uii" are given by 



Li(a) = Xi-i(a), 

Si(a) = -(Xi{a) ~ Xi-i(a)), 

Xi{a) = -giXi(a), 

Xi{a) 



(5.26a) 
(5.26b) 
(5.26c) 
(5.26d) 

(5.26e) 



Since the resonator and qubit frequencies are pulled by 
the classical field due respectively to the nonlinearity and 
the ac-Stark shift, the Lamb shift depends on these pulled 
frequencies, and therefore on the amplitude of the cavity 
field. 



Finally, projecting the qubit onto its {|0) , |1)} subspace 
and tracing out the resonator degrees of freedom yields a 
reduced qubit master equation 

i" 

Pq = ~i [H,p q ]+7±' V[a-]pq+7f' V[cr+]p q + -7j-'D[(T z ]p q . 

(5.27) 

In this expression, we have defined 

bj'" 

H = -fa z + .g (ao, s e-^V_ + h.c), (5.28) 



where uj"q = u)'{'(a) — uj'^'(a), and 



where 



Tip 7ip "b r^m, 



_ nD 2 knl|c*i— «ol 



Lpm 



^* 7|2X y d - Xfa 14 \ 2 



(5.29) 



(5.30) 



and D = \ai — a \ is the distance between the pointer 
states. In the equation for the effective dephasing rate 
r vm , we see the measurement-induced dephasings due to 
single-photon cavity losses (second term), to two-photon 
cavity losses (third term), and to the information car- 
ried out by the excitation emitted when the qubit re- 
laxes. While these three channels leak information about 
the qubit state, only the single-photon cavity loss chan- 
nel is usually monitored. Moreover, since in practice 
k >• -fX 2 /g 2 , knl, only this last channel will convey any 
significant amount of information and contribute to qubit 
dephasing. 



VI. BACKACTION ON THE QUBIT 

Following the reduced qubit model derived in section [Vj 
here we revisit the results presented in section III for the 
dispersive regime of linear circuit QED. In this section, we 
compare the theoretical model to experimental data and 
numerical simulations. The parameters used throughout 
are given in the caption of Fig. [4j These parameters 
were adjusted to fit independent spectroscopic and time 
domain measurements of the device used in Ref . [16] . This 
device was composed of a transmon qubit [3] coupled 
to a coplanar waveguide resonator made nonlinear by a 
Josephson junction embedded in its central conductor. 

In subsection |VI A| we first quickly present the experi- 
ment already described in Ref. [TS]. We then look more 
precisely at the Lamb and ac-Stark shifts of the qubit 



transition frequency w 10 in subsection VI B and at its 
linewidth in subsection IVI CI 



A. Experiment and qubit spectra 

In Ref. [16], we presented spectroscopic measurements 
of a transmon qubit coupled to a driven nonlinear res- 
onator. The qubit was probed through the resonator 
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with a drive of amplitude e s and frequency u> s ~ uiifi. 
Meanwhile, that resonator was pumped with a drive of 
amplitude e p and frequency uj p ~ w,.. The pump field 
was applied long before the qubit probe was turned on, 
enabling the resonator to reach its stationary state. Two 
detunings between the pump frequency lo p and the res- 
onator frequency uj r where studied in detail. This was 
done in order to explore both the parametric amplifi- 
cation and the bifurcation regimes. Consequently, two 
biasing points ui p /2n = (6430,6450) MHz corresponding 
to Vt/tt c — (3.1,0.7), are presented below. Here, we rede- 
fined Q with respect to the effective resonator frequency 
as pulled by the qubit in the ground state rather than 
the bare resonator frequency. These two biasing points 
are illustrated by the two vertical lines in the stability 
diagram of Fig. [3] (b). After probing the qubit, a bifurca- 
tion measurement was performed in order to determine 
the probability -P(|l)) that the qubit was excited by the 
probe drive. 

The resulting experimental spectra are presented in 
the top panels of Fig. El for Q/Qc = 3.1 (top left) and 
Q/Qc = 0.7 (top right) as a function of the pump drive 
amplitude. The pump amplitude (horizontal axis) is 
converted to a logarithmic scale to match the experimen- 
tal power in decibels, up to a constant offset that was 
calibrated in Ref. [16 . In the bifurcation regime (top 
left, n/Qc = 3.1), we clearly see the jump in the qubit 
frequency associated with the jump from the low am- 
plitude to the high amplitude dynamical states of the 
resonator. We also see that the line remains narrow and 
actually tends to narrow down at higher powers. In the 
parametric amplification regime (top right), we see a 
more monotonous shift of the qubit line with the mea- 
surement power with an important broadening around 
201og(e p /27r) = 22. 

These spectra are then compared to the analytical 
steady-state solution of the reduced qubit master equa- 
tion Eq. (5.27) in the bottom panels. The exact analytical 



Experiments 



solution of this equation yields 



P(|l» = 



7^(7 2 2 + a 2 )+ 2 7 2|goao, s | 2 
(7f: o +7r:o)(7 2 2 + 5 2 )+4 72 | 5o ao,. 



where 



72 = V" + 



7j.,o 



7f,o 



(6.1) 



(6.2) 



; i,o 



When comparing the experimental to the analytical 
spectra, we notice small deviations between the back- 
ground level as well as the amplitude of the spectroscopy 
lines. Aside from the limits of our model, three effects 
can cause these deviations. First, there is experimental 
thermal noise — which should not exceed 50 mK — that 
is not taken into account in the theory and may yield a 
minor thermal qubit excited state population. Second, 
the experimental excited state population is extracted 
from the probability of bifurcation, which can yield an 
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201oi 



(e p [MHz]/2^) 



FIG. 4. (Color online) Experimental (c.f. Ref. [16]) 
and analytical qubit excited state |1) population 
for the two operating points indicated in Fig. [3] 
The qubit is a transmon with bare parameters 
(w 1>0) W2,i, 7,7 v )/27r = (5720,5421.6,0.22,0.25) MHz. 
The resonator's bare parameters are (uj r , K, K' , k, knl)/27t = 
(6453.5, -0.625, -0.00125, 9.6, 0) MHz and the qubit-resonator 
couplings are (go,gi)/2-K = (42.4,58.4) MHz. These param- 
eters were chosen to fit those of Ref. [T6 1. Couplings to 
higher transitions as well as higher transition frequencies 
can be computed from the transmon Hamiltonian [6]. The 
experimental attenuation required to link the experimental 
power in dB to the theoretical parameter e p was calibrated 
in Ref. |16j . Top: experimental spectroscopy results from 
Ref. [16]. Bottom: analytical stationary solution Eq. (6.1 1. 
Left (Right) : pump frequency w p /2tt = 6430(6450) MHz, 
corresponding to 0,/Q.c = 3.1(0.7). These points are identified 
by vertical lines on Fig. [3] The amplitude £ s /2tt = 3 MHz 
was chosen to fit the experimental power broadening of the 
qubit lines. 



error of at most 0.05 in the estimated population. Third, 
the correspondance between the theoretical amplitude e s 
of the spectroscopy drive and the experimental amplitude 
could not be calibrated as precisely as the calibration 
provided by the ac-Stark shift for the pump drive [To] . 
Overlooking these deviations, other experimental features 
such as the spectroscopy lines' position and width are 
qualitatively reproduced by our analytical spectrum. In 
the following sections, we quantitatively compare these 
to our model. 



B. Lamb and ac-Stark shifted qubit frequency 

The experimental spectra presented in Fig. [4] were fit- 
ted using Lorentzian and the peak positions and widths 
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FIG. 5. (Color online) Lamb and ac-Stark shifted qubit fre- 
quency as a function of the drive strength e p for the two operat- 
ing points indicated in Fig.[3]and used in Ref. [16]. Parameters 
are the same as Fig|4] Points are experimentally (black circles) 
and numerically (orange squares) extracted qubit transition 
frequency uji.o . Lines are analytically computed w"' — oj'q 
with the complete Eq. (5.26a l (full black lines), when setting 

= (dotted red lines) or when taking lo p = lo t in S? and 

(dashed green lines). 



were extracted from those fits, yielding the qubit transi- 
tion frequency and dephasing rate. We also numerically 
integrated the multi-level Jaynes-Cummings master equa- 
tion (2.6) to obtain numerical spectra that were fitted 



using the same procedure. The qubit frequency extracted 
from experimental (black circles) and numerical (orange 
squares) spectra is plotted in Fig. [5] as a function of the 
pump power for the two operating points. Numerical 
simulations and experimental data almost coincide, sug- 
gesting that the initial master equation (2.6 1 contains all 
the relevant physics. 

We then compare these data points to three versions 
of the dispersive approximation. Full black lines corre- 
spond to the complete equation (5.26a I, dotted red lines 



correspond to the second order approximation for the 
dispersive shift (i.e. = 0), and dashed green lines cor- 
respond to setting uj p = uj r when calculating §?■ and K?. 
Since the parametric amplification regime (right panel) 
correspond to a pump drive very slightly detuned from 
the resonator frequency, as well as to a low number of 
photon (n ~ 20), all three curves almost coincide in this 
regime. 

On the other hand, in the bifurcation regime (left panel), 
both the pump-resonator detuning and the number of 
photons after bifurcation are larger (n ~ 50), yielding 
a significant difference between the three curves above 
bifurcation. We see that the assumpti on lo p — u r (dashed 
green lines), which as discussed in Sec. IV A| is often made 



when calculating the ac-Stark shifts, yields a shift that 
is too small. This is expected since assuming lj p = u) r 
yields a larger qubit-pump detuning, and correspondingly 
smaller values of Sf and Kf. This effect can also be 
confirmed at lower power although it is not visible in these 
plots. We also see that the second order approximation 
(dotted red lines) yields a dispersive shift that is too 



large. This is also expected since the sign of each order 
in perturbation theory alternates sign in the dispersive 
regime and since the fourth order is contained in the full 
model. 

With this model, the qubit can be used as a tool to 
characterize the nonlinear resonator. Indeed, the distance 
between the resonator's low and high amplitude states 
at the threshold of bifurcation directly depends on the 
resonator nonlinearity K and the drive frequency u> p and 
amplitude e p . While experimentally u p is known to a very 
high precision, the resonator nonlinearity K can only be 
estimated to about ±30% from the design parameters due 
to its nonlinear dependence on sample parameters [16] . 
Moreover, the experimental line attenuation A between 
the source and the input of the sample — which is required 
to make the correspondance between the experimental 
power P p and the theoretical parameter e p — can only be 
estimated up to about 2 dB [16] . Performing a series of 
spectroscopic measurements for many pump frequencies 
ujp and fitting the extracted qubit frequencies to the model 
derived here then makes it possible to extract both K and 
A with improved precision. This was done in Ref. |16j 
and resulted in an uncertainty of 2.4% for K and 0.2 dB 
for A] a ten-fold improvement in precision. 



C. Qubit linewidth and validity of linear response 

We now examine the linewidth of the qubit transition. 
We know that, in addition to the intrinsic dephasing rate 
72,int = 7(p+7/2, the lines are broadened by measurement- 
induced dephasing [25] and by dressed-dephasing g3]- I n 
addition, there is always some power broadening due 
to the finite spectroscopy power. Here, we are mostly 
interested in the measurement-induced dephasing and how 
it is modified by the nonlinear nature of the resonator. 
The experiments presented in Ref. [TB] and whose results 
are reproduced here were therefore carried in a regime 
where power broadening is small. Moreover, since there is 
no dependence of the experimental background population 
over the pump power, we assume that dressed-dephasing 
is also negligible due to a small amplitude of dephasing 
noise at GHz frequencies. The only additional dephasing 
source is therefore measurement-induced dephasing in 7"' 



given in Eq. ( 5.29 ) and in practice is dominated by the 



k\ol\ — ao\ /2 contribution. 

We present in Fig. [6] the half- width at half-maximum 
of the spectroscopy lines as a function of the pump power 
for the two operating points fi/Qc = 3.1(a), 0.7(6). Grey 
circles (orange squares) are again the widths extracted 
from experimental (numerical) data. Full black lines are 
the analytical widths 72/2^ given by Eq. (6.2). Dashed 
green lines are the same as the full black lines, but using 
linear response theory for the fields on %p instead of the 



solutions of Eq. (5.19). More precisely, we obtained the 



dashed green lines taking 



if) + 3K\a\ 



(6.3) 



12 



10 



(a) ' 1 h/n c = 3.1 1 


(b)i ' ' n/b c = 0.7 1 . 
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FIG. 6. (Color online) (a) and (b) : Qubit line's half-width at 
half-maximum as a function of the drive strength e p for the 
two operating points indicated in Fig. [3] and used in Ref. [16] . 
Points are experimentally (grey circles) and numerically (or- 
ange squares) extracted qubit linewidths. Lines are analytical 
solutions corresponding to 72 where the fields Uj are com- 
puted according to the nonlinear response Eq. (5.191 (full 



black lines), to linear response Eq. (6.3| (dashed green lines) 
and by replacing k|qi — Qo| 2 /2 with the linear resonator result 
Eq. (6.4| (dotted red lines), (c) and (d) : Phase space repre- 
sentation of the fields ao (black lines, circles) and on (red lines, 
squares) as given by the linear (dashed lines, empty symbols) 
and nonlinear (full lines, full symbols) response theories. All 
four symbols of a given set (1 or 2) correspond to a given 
pump amplitude e p . 



= 0. 



where a is the solution of Eq. (5.19) with = 
Finally, dotted red lines we obtained by replacing T vm 
by the result of Ref. 28J for a linear resonator 



k 2(\ ai J 2 + K,p| 2 )x 2 

2 K 2 /4 + X 2 + (W r - U p ) 



(6.4) 



■Sg. 



where X = §1 

The first striking observation is that, contrary to circuit 
QED with a linear resonator [32]) the linewidth does not 
strictly increase with the drive power or equivalently 
with the number of photons in the resonator. In fact, in 
the bifurcation regime [Fig. [6]Ja)], the linewidth shows 
a sharp maximum at the bifurcation power, whereas in 
the parametric amplification regime, the linewidth shows 
a smooth maximum at a power that corresponds to the 
maximum gain of the amplifier |16j . This is illustrated 
by the lack of even qualitative agreement between both 
experimental and numerical data points and the result 
expected for a linear resonator (dotted red line). 

Narrowing of the linewidth at high power is predicted 
both by the nonlinear (full dark lines) and the linear 



(dashed green lines) response theory. However, while both 
give a qualitative agreement with experimental and numer- 
ical data points, only the nonlinear response theory gives a 
quantitative one. In the bifurcation regime (Q/fic = 3.1), 
the nonlinear response theory reproduces the experimen- 
tal behavior with good accuracy on the whole range of 
powers, whereas linear response predicts bifurcation at too 
low power and linewidths twice as large at bifurcation. In 
the parametric amplification regime (fi/fic = 0.7), only 
the nonlinear response solution gives semi-quantitative 
agreement near the maximum linewidth, while linear re- 
sponse theory predicts a much lower linewidth. How- 
ever, even the nonlinear response solution mispredicts the 
linewidth when it is above ~ 5 MHz. We explain this by 
the breakdown of the \ai — oiq\ < 1 approximation, which 
corresponds to a measurement-induced dephasing rate of 
about T vm ~ k/Att ~ 5 MHz. 

To understand the non-monotonous behavior of the 
linewidth with drive power, we refer to Figs. [6] (c) and 
(d), where we plot the value of the fields &o(i),p as black 
(red) lines in the complex plane for the two operating 
points, for a range of power e p /2ir £ [0, 150] MHz and for 
nonlinear (full lines) and linear (dashed lines) response 
solutions. We see with these plots that even though the 
number of photons increases as the distance to the origin 
grows, the distance between the solutions a% tP and arj,p 
does not. In fact, the distance D can be as small at higher 
power than at small power. 

For reference purposes, we also plot two sets of four 
points in pannels (c) and (d). Each set corresponds to a 
given pump amplitude e p , for nonlinear (full symbols) and 
linear (empty symbols) theory, and for ao (black circles) 
and o?i (red squares). Comparing the points within a 
given set of four points, we can see that a larger distance 
between a circle and its corresponding square — and 
hence the larger the gain of the amplifier — correspond 
to a larger disagreement between the linear and nonlinear 
solutions (distance between a full and a corresponding 
empty symbol). 

We can compute a range of validity of the linear re- 
sponse theory by computing the fields a, lP to second 
order (i.e. quadratic response theory). If we define 



.,(!) 

i .p 



* (2) 



where <x- 2 is the second term of 

„(2 



equation Eq. ( 6.3 ) and a\ is the next order correction, lin- 



ear response theory will be valid if the ratio r = af^/af^ 
is small. Since for a qubit measurement, the signal that is 
amplified is a frequency shift S = ±(§i — So), we can de- 
fine a maximal value of S that allows r to be smaller than 
a threshold r t in the region of highest gain. This maximal 
value 5 mal , computed using a conservative value of 10% 
for the ratio of the quadratic correction over the linear 
correction, is plotted in Fig. [7] as a function of the reduced 
detuning Q/Clc- We see that the maximal coupling for 
the parameters given in the caption of Fig. [4j typical for 
circuit QED, never exceeds about 0.5 MHz. Moreover, 
the maximal coupling in fact vanishes when approaching 
the critical detuning Qc- This maximal coupling is to 



13 



1000 



N 
CN 



CO 




0.5 



FIG. 7. Maximum dispersive coupling S max for which the 
linear response theory is valid with 10% uncertainty at the 
region of maximum gain as a function of the reduced detuning 
Sl/Slc in the parametric amplification regime. 



be compared with the resonator linewidth k in order to 
determine if it is viable for a qubit measurement. With a 
realistic criteria of x > 0.2k to get a good measurement, 
one therefore needs either k/2tt ~ 1 MHz or a smaller 
nonlinearity K in order for linear response theory to be 
valid in this system. The former however implies a longer 
measurement time, while the latter implies a smaller gain, 
both impairing the efficiency of the measurement. It 
therefore seems unlikely that linear response theory will 
be sufficient to describe any superconducting qubit read- 
out using a nonlinear resonator until the qubit lifetimes 
become long enough for longer measurement time to be 
viable. 



D. Quantum limit to the added noise 

Using the results presented in Fig. [6| we can try to 
answer the question of whether or not a dispersive homo- 
dyne measurement using a nonlinear resonator can reach 
the quantum limit T vm — T meas /2 as is the case for a 
linear resonator [27] . Indeed, assuming small squeezing, 
if one were to make a homodyne measurement using the 
pump drive, the measurement rate would be given by 
= k\(Xi — «o| 2 HZ]- Since this measurement rate 
is exactly twice the dominant part of the measurement 
induced dephasing caused by these same pump photons 



r„ m given at Eq. (5.301, we can say that the quantum 



limit is reached if the theoretical prediction fits the exper- 
imental linewidth. If the experimental linewidth is larger 
than the theoretical prediction, it however means that the 
limit is missed. Finally, if the experimental linewidth is 
smaller than that predicted by the model, it means that 
one of the approximation is probably breaking down. 

Looking at Fig. [6] (b), we then reach a different con- 
clusion whether we consider linear or nonlinear response. 
Indeed, around 20 log(e p /27r) € [20,30], the experimental 
linewidth is much higher than the prediction from lin- 



ear response, and we would therefore conclude that the 
quantum limit is missed by the measurement. This is 
qualitatively the same conclusion as the one obtained by 
Laflamme and Clerk [55], also in a linear response theory. 
However, we know from Fig[7j that for Sl/Clc — 0.7 as in 
Fig.[6](b), the maximum dispersive coupling supported by 
a linear response treatment is S max /2n ~ 200 kHz, about 
four times smaller than the one used here. If we now com- 
pare the nonlinear response model prediction (black line) , 
we see that it matches the experimental observations on 
a much wider range, and we recover the quantum limit in 
this range. There is also a regime where the theoretical 
prediction is above the experimental observation. This 
regime corresponds to a linewidth ~ kD 2 /2 > 5 MHz 
since k/2tt ~ 10 MHz, and therefore to D > 1, breaking 
the small distinguishability approximation that we have 
made. Therefore, while our result shows that the quantum 
limit can be reached with a nonlinear resonator, the ques- 
tion remains open in the case of large distinguishability 
or large squeezing where our model breaks down. 



VII. CONCLUSION 

In summary, we have derived an analytical model 
to describe the backaction of a driven nonlinear res- 
onator on a multi- level qubit. This is done using unitary 
transformations, and especially using the polaron |25F 
[27] and dispersive [23] [24] transformations. We obtain 
a reduced model that contains the physics of the lin- 
ear and quadratic ac-Stark shifts as well as the Lamb 
shift of the qubit frequencies. The model also contains 
dressed-dephasing [23j |43l 04] , Purcell relaxation [8] and 
measurement-induced dephasing [571 HH1 B2] ■ Contrary 
to other theoretical models, both qualitative and quanti- 
tative agreements are found for the ac-Stark and Lamb 
shifted qubit transition frequencies as well as for the qubit 
linewidth. 

Moreover, the model that we have derived here goes 
beyond some assumptions that are frequently made and 
that are valid in the case of a driven linear resonator, 
but not in the nonlinear case. These assumptions are the 
resonant driving of the resonator, the linear response of 
the resonator field to the qubit signal and the two-level 
character of the qubit. Considering detuned driving of the 
resonator yields linear and quadratic ac-Stark shifts that 
depend on the qubit-drive frequency detuning rather than 
the qubit-resonator frequency detuning and are therefore 
slightly different than usual dispersive shifts [5] . Going be- 
yond linear response theory yields measurement-induced 
dephasing rates that are qualitatively different from those 
found with linear response and that are found to match 
the experimental and numerical data in most regimes 
considered. In particular, we show that the measurement- 
induced dephasing rate does not increase with the mea- 
surement power or the number of photons, but rather 
with the distance between two pointer states ot\ and ao of 
the resonator fields. The precise quantitative agreement 
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between the model and the experiment has also allowed 
us in Ref . [16j to characterize the nonlinearity of the res- 
onator and the attenuation of the transmission line with 
an accuracy ten times better than what was otherwise 
achievable. 

We have finally also shown that the results given by 
linear reponse theory are unlikely to apply to any high- 
fidelity qubit measurement using a nonlinear resonator. 
One consequence of this is to reopen the question of 
whether or not measurement with a nonlinear resonator 
is quantum limited in the amount of dephasing it causes 
on a qubit. Indeed, while Laflamme and Clerk [53] have 
shown that the quantum limit is missed by a factor G, 
the gain of the amplifier, this result was obtained in a 
linear response theory and therefore is not applicable in 
the systems considered here. This question then remains 
open and could be answered using a quantum trajectory 
approach as was done before for a linear resonator |27j . 
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Appendix A: Polaron transformation 

In this appendix, we give the result of applying the 
polaron transformation (5.3 1 on the different parts of the 



master equation (2.6 1. The building blocks from which 



all operators can be transformed are 

a' = a + Tl a , 
n^j+i = n M+ iD t (a l )£>(a. i+ i) « H M+ i, 
TT- ■ — TT 



(Ala) 
(Alb) 
(Ale) 



where we noted the transformed-frame operator O' = 
P^OP. Using these relations, transforming the Hamilto- 
nian H r yields 



H' r = w r {a)a + U* a a + H a aJ + \U a \ 2 
K' 



K 



jn a | 4 + (2a t |n Q | 2 n Q + h.c.) + Aa)a\Ii a \ 2 -I- (a^U 2 a + h.c.) 



n Q | 6 + (3at|n Q | 4 n a + h.c.) + g^air^i 4 + (3^ \u a \ 2 u 2 a + h.c.) 



(A2) 



where we have dropped terms with more than two res- 
onator ladder operators. This approximation assumes 
that \a\ 3> I (a) I in the transformed frame. We will see 
that with a proper choice of a, the resonator in the trans- 
formed frame is close to its ground state. Transforming 
H q is trivial since it is diagonal in the qubit subspace 
and therefore commutes with the transformation and 
H' = H q . Transforming the interaction Hamiltonian Hj 
yields 

M-2 

H'j=Y^9i W + «] [ n «+i + Hi+i,<] 



i=0 
M-2 



(A3) 



E 9* [n* it ,t+i + iF+i^ria] . 



i=0 



In obtaining this equation, we assumed that D{on+\ — 
aj) « 1 and made a RWA for the second line. Not doing 
the RWA would only yield a small Bloch-Siegert shift to 
the qubit transition frequencies [S5]. However, we choose 
not to do a RWA on the first line at this point. This 
will allow us to get the sidebands Hamiltonian for a MLS, 
equivalently to what was done in Ref. [3H] for a TLS. 
Transforming the drive Hamiltonians H s and Hd is also 



trivial and yields 

H' d = E ^"""V + H^) + h.c, (A4a) 

d 

H' s = e s e- iw *V + H*) + h.c. (A4b) 

Finally, since the transformation P moves the system to 
a time-dependent frame, a Hamiltonian 

Hp = iP^P = (-TLV + m* a a) - Ina[n a n*], (A5a) 

must be added to the total Hamiltonian in the transformed 
frame. 

For the dissipation, we can show that 



V[a']p' = V[a]p' 



-n> + h.c,p' 



(A6) 



y,n*] + h.c. 



In this equation, the second term is of Hamiltonian form 
and will be added to the Hamiltonian in the transformed 
frame. It is worth noting that, if p' is the ground state of 
the resonator in this frame, the last line is equal to zero. 
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For the two-photon dissipation, we get 



v[ a ' a y « D[n 2 >' + 4p[ a n a ] P ' 
-» ^a|n Q | 2 n; + i 
+ (2an Q + a 2 ) [p',n; 2 ] +h.c, 



i Ln Q | 2 n; + in; 2 a 2 ') + h.c. )/9 ' 



(A7) 



where again, the last line is zero if p' is the ground state of 
the system, and the second line is of Hamiltonian form and 
will be included in the Hamiltonian in this transformed 
frame. When obtaining this result, we assumed again 



that | a | 3> | (a) |, and neglected any term with more than 
two ladder operators. 

Since the polaron transformation commutes with II e , 
2?[II £ ] stays the same in the transformed frame. Moreover, 
since we assumed that |atj+i— a<| < 1, we do not transform 
the dissipators DpI^i+i]. 



Appendix B: Dispersive transformations 

In this appendix, we transform the different parts of the 
master equation ( 5.10 1 according to the classical dispersive 



transformation Dp given at Eq. (5.12). We note O 



D^O'Dc*. To first order in £, we get 



M-2 



M-2 



i=0 ' i=0 



(Bl) 



where 



M-2 



iJ Dc = it> c T>l « i Yl ~ 'til,. (B2) 



Assuming the form of Eq. ( 5.13 ) for we can compute 
£,(t). Doing this and taking 

J2* K w <+M ~ u d )^,d + g l a % , d ] e~ luJdt = 0, (B3) 



or equivalently 



-giott.d 



makes the non-diagonal terms in Eq. ( |B1[ ) vanish. With 
this choice, transforming H' to fourth order in perturba- 
tion theory and assuming large frequency differences and 
sums \uJdi ± uJd 2 \ as well as small \oti t d — Oi+i,d\ yiel ds the 
Hamiltonian Htf = T^ c H'qD c + # Dc given by Eq. (15 . 14h . 



In order to transform H[, we need to know how to 
transform a diagonal operator 



M-2 



D^n,D c w il x + ( x t+i ~ XiMUi+i,i + e? n< (i +i), 

(B5) 



i=0 



(B4) and the off-diagonal operator £_ = J2i=o 2 .9;IT,i+ 



M-2 



d^e d c » 5» n M+i + ns*n ai e- toit + i 5] n Kd |n ad | 2 n Qd e- 



iuidt 



i=0 



M-2 



i E E A * 



M-2 



-^i+^^n^i + h. 



i £ £ [ 2A * 1 + X t 2 X ti) - ft (Af-i^ + A&A* x 

i__0 di.da 



a^a^e^-^^+x+h.c. 

(B6) 



where we made the same assumptions as previously and Xf, Af, Sf and are defined at equations (5.16)-(5.17) 



These two equations can be combined and used to transform H[ yielding the result of Eq. (5.18) with 

G" = C + fn sd + ^n Kd |n Q | 2 ) n ad < 



-iu d t 



(B7) 
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and 



M-2 



2=0 di,d2 
M-2 



a dl a d2 e~^+^n i+hl + h.c.] } | 



E E { [ 2A * + X t 2 - X ti) ~ 9, (A-ti^ti + A m A m)] [a* di a d2 e^-^n^ +1 + h.c.] } \ 

i=0 dijda 
M-2 

E A M^i -§- 2 ] [« dl a d2 e- l ^i+^)*n i+1 , + «* i a d2 e 4 ^i^)*n i , +1 ] t 
+ h.c. 

(B8) 



In this last Hamiltonian, the choice of the polaron frame Eq. (5.191 has already been made. Finally, transforming 
dissipators according to D c yields 



M-2 



2 



ti 



p' -> 2 7v r> 



5i 



M-2 



2 ^ E E* * 2 \^i\\mi+iAp" + ^[Hm+i]/) (B9a) 



i=0 d 



9o 



M-2 , x 2 

9i 



M-l 

E 



9o 



AI — 2 



kV [n Q ] p' -> K p [n a ] P " + K ^ 5^ |a i+1 - ai\ 2 \^ dti \ 2 (p[n i+1 ^p" + V[U hl+1 }p"). 



(B9b) 



(B9c) 



i=0 d 



These transformed Hamiltonians and dissipators can be combined to obtain the master equation in the polaron and 
classical dispersive frame given in Eq. (5.21 1. 
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